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Abstract. The idea to use classical hypergeometric series and, in particular, well- 
poised hypergeometric series in diophantine problems of the values of the polylogarithms 
has led to several novelties in number theory and neighbouring areas of mathematics. 
Here we present a systematic approach to derive second-order polynomial recursions 
for approximations to some values of the Lerch zeta function, depending on the fixed 
(but not necessarily real) parameter a satisfying the condition Re (a) < 1. Substitut- 
ing a = into the resulting recurrence equations produces the famous recursions for 
rational approximations to C(2), C(3) due to Apery, as well as the known recursion for 
rational approximations to £(4). Multiple integral representations for solutions of the 
constructed recurrences are also given. 



The idea to use classical hypergeometric series [9, 11] and, in particular, well-poised 
hypergeometric series [13] in diophantine problems of the values of the polylogarithms 
has led to several novelties in number theory and neighbouring fields of mathematics. 
Here we present a systematic approach to derive second-order polynomial recursions 
for approximations to the numbers 
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where the fixed (but not necessarily real) parameter a satisfies the condition Re(a) < 1. 
Substituting a = into the resulting recurrence equations produces the famous recur- 
sions for rational approximations to £(2) = 2Z^"(0), £(3) = Z 3 (0) due to Apery [1], 
as well as the recursion for rational approximations to £(4) = Z^O) known as the 
Cohen-Rhin-Sorokin-Zudilin recursion ([8, 15, 17]), which is proper from both the 
historical and the alphabetic point of view. 

To make clear to the reader, what do we mean by a recursion for approximations 
to a number z £ C, we introduce a formal definition. The requirement to such the 
recursion is to have two linearly independent solutions {u n }'^ =0 and {f n }^L (uniquely 
determined by the recursion itself and initial conditions) such that v n /u n —> z as 
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n — > oo. In the case z e K (e.g., corresponding to a G Q in the above definitions), 
we usually restrict ourselves to the sequences {«„}™ =0 , {v n }^ =0 consisting of rational 
numbers only and may interpret them as denominators and numerators, respectively, 
of rational convergents to z. 

We apologise in advance for facing the reader with, sometimes, cumbersome for- 
mulae. Although the ideas of the well-poised hypergeometric construction of linear 
forms r n = u n z — v n , n = 0, 1, 2, . . . , are simple, the appearance of lengthy formulae is 
unavoidable in this type of analysis. Our main results are the recursions (4), (8), (10), 
(12), (13), and the integral representations (6), (9), (11) (generalizing those of [3]) 
for the linear forms {r n }^ =0 . 

1. Z^ipt). Take the rational function 

fi„( i ) = n !.(2 i -2a. + n )(^ 1) ''' (i -" ) ' (i - 2Q + " +1) ''' (i - 2a + 2 " ) 
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T(t-a + n+l)J T(t-2a + n + l) 
satisfying the property 

t^-(t-2a + n): R n (t) i-> (-l) n R n (t). (1) 
Decompose the function R n (t) as the sum of partial fractions, 

j=0 fc=0 v ' 

and use the last representation as in [19], proof of Lemma 1, to sum the quantity 
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The property (1) yields UQ n = U2 n = for all n = 0, 1, 2, ... , hence setting it n = u\ n j2 
we obtain the linear forms 

r n = u n -2Z^(a)-v n eQZ^(a)+Q, n = 0,1,2,..., (3) 
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with effectively determined coefficients u n and v n . Applying Zeilberger's creative 
telescoping [12] in the manner of [19], Section 2 (namely, computing the certificate 
and the correponding difference annihilating operator for the function (— l)*i? n (t), 
which is rational with respect to either t or n), we arrive at the following recursion 
satisfied by both the linear forms (3) and their coefficients: 

(n + l) 2 (n + 1 — a) 2 (5n 2 — Aan + a 2 )-u n _|_i 

- (55n 6 - ll(14a - 15)n 5 + (179a 2 - 385a + 180)n 4 

- (116a 3 - 358a 2 + 332a - 85)n 3 + (45a 4 - 174a 3 + 232a 2 - 113a + 15)n 2 

- a(a - l)(10a 3 - 35a 2 + 41a - 12)n + a 2 (a - l) 2 (a 2 - 3a + 3))u n 

- n 2 (n - a) 2 (5(n + l) 2 - 4a(n + 1) + a 2 )w n _i = 0, (4) 

uq = 1, u\ = a — 3a + 3, 

a 2 -4a + 5 lim — = 2Z 2 "(a). 

Vq =0, 17! = — n ^oo Un 

(a — ±y 

(We justify the latter limit relation by notifying the following consequence of the 
forthcoming formula (6): r n — > as n — > oo.) 

On the other hand, the series (2) may be easily identified with the very-well-poised 
hypergeometric gF 5 (— l)-series 



r n = (-l) n nl J2( 3n + 2 - 2a + 2z/) 

r(3n + 2 - 2a + u) Tjn + 1 - a + u) 3 T(n + 1 + u) 
X r(l + z/)r(z/ + n + l-a) 3 r(2n + 2-2a + z/) ^ ' ' ^' 

which admits the double integral representation 

r n = u n - 2Z 2 -(a) - v n = (-1)" / / (1 _ x(1 - dxdy ^ 
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(see [17], Theorem 5). Whipple's transformation ([2], Section 4.4, formula (2)) gives 
one a more direct way to deduce the integral (6): first convert the series (5) into the 
hypergeometric 3_F2(l)-series 

n [U+ a) ^ Q r(l + z/)r(2n + 2-a + z/)2 

and secondly use the Euler-type integral formula for the latter series. 
2. Z 3 (a). This time, take the rational function 

R n (t) = • (2t - 2a + „ ) WljN^±iL±^ 

((* - a)(t - a + 1) • • • (t - a + n)) 
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satisfying the property 

t ■->-(*- 2a + n) : R n (t) i-> -R n (t). (7) 
After partial-fraction decomposition we arrive at the quantity 

- R n(t) = u n Z 3 (a) -v n E QZ 3 (a) + Q, n = 0, 1, 2, . . . , 



t=i 



with effectively computable coefficients u n and v n . Zeilberger's creative telescoping 
produces the recursion 

(n + l) 3 (n + 1 - a) 3 (2n - a)(3n 2 - 3an + a 2 )w n+ i 
- (2n + 1 - a) (I02n 8 - 408(a - l)n 7 + 2(359a 2 - 714a + 321)n 6 

- 6(a - l)(121a 2 - 238a + 83)n 5 + 3(152a 4 - 605a 3 + 811a 2 - 415a + 64)n 4 

- 2(a - l)(89a 4 - 367a 3 + 461a 2 - 177a + 15)n 3 
+ a(40a 5 - 267a 4 + 634a 3 - 669a 2 + 304a - 45)n 2 

- a 2 (a - l)(2a - l)(2a 3 - 17a 2 + 37a - 25)n - a 3 (a - l) 2 (2a 2 - 6a + 5))u n 
+ n 3 (n - a) 3 (2(n + 1) - a)(3(n + l) 2 - 3a(n + 1) + a 2 )w n _i = 0, (8) 

uq = 1, u\ = 2a 2 — 6a + 5, 

v 

(a 2 - 3a + 3)(a - 2) hm — = Z 3 (a), 

V = 0, Vl = -. n ^°° U n 

(a — l) 6 

while writing r n as a very- well-poised hypergeometric s-FY(l)-series and applying [17], 
Theorem 5, we obtain the triple integral 

r n = u n Z 3 (a) - Vn = l - [[[ ' ^"'f " ( ' d* dy d, (9) 

2 77,/ (1 - x(l - y(l - z)))^ 1 y y > 

[0,1] 3 

Bailey's transformation ([2], Section 6.3, formula (2)) allows us to write the quan- 
tity r n as the Barnes-type integral 

_ f c+lo ° T(n + 1 + s) 2 T(n + 1 - a + s) 2 r(-s) 2 d ^ 
Tn ~ Ani J c _ loo T{2n + 2 - a + s) 2 S 

1 ^ d ( (t- l)(t- 2) •••(£- n) 



2 f-^dt \(t-a)(t-a + l)---(t-a + n) 



where the real constant c lies in the interval Re(a) — 1 < c < 0. This gives another 
way to deduce the recursion (8), by applying Zeilberger's creative telescoping directly 
to the latter summation (cf. [18], Lemmas 1-3, for the proof in the particular case 
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3. Z±(a). Finally, take the rational function 




satisfying the property (7) and consider the quantity 




-^r~ = Un ' 6Za ^ ~ VnE Q ^ 4(a) + Q ' = 0, 1, 2, . . . . 



Then Zeilberger's creative telescoping gives the recursion 



(n + lf(n + 1 - a) 3 (n + 1 - 2a)(39n 4 - 65cm 3 + 45a 2 n 2 - 15a 3 n + 2a 4 )u n+1 

- (l0530n 13 - 1755(40a - 39)n 12 + 18(11881a 2 - 23400a + 10881)n n 

- 9(43964a 3 - 130691a 2 + 122499a - 36036)n 10 

+ (497482a 4 - 1978380a 3 + 2795153a 2 - 1651455a + 343161)n 9 

- (449452a 5 - 2238669a 4 + 4229444a 3 - 3756546a 2 + 1559025a - 241137)n 8 
+ 2(149999a 6 - 898904a 5 + 2128142a 4 - 2523748a 3 

+ 1567577a 2 - 480285a + 56394)n 7 

- (149336a 7 - 1049993a 6 + 2995163a 5 - 4449872a 4 
+ 3679649a 3 - 1676024a 2 + 385125a - 33930)n 6 

+ (55088a 8 - 448008a 7 + 1503025a 6 - 2693161a 5 + 2786514a 4 

- 1681907a 3 + 568968a 2 - 96291a + 5967)n 5 

- (14696a 9 - 137720a 8 + 536294a 7 - 1132580a 6 + 1413762a 5 

- 1065166a 4 + 474344a 3 - 116539a 2 + 13455a - 468)n 4 

+ a(a - 1) (2692a 8 - 26700a 7 + 105832a 6 - 220076a 5 + 260191a 4 

- 176174a 3 + 65540a 2 - 11955a + 780)n 3 

- a 2 (a - l) 2 (304a 7 - 3430a 6 + 14198a 5 - 29252a 4 
+ 32370a 3 - 18825a 2 + 5265a - 540)n 2 

+ 2a 3 (a - l) 3 (8a 6 - 128a 5 + 581a 4 - 1198a 3 + 1220a 2 - 558a + 90)n 
+ 4a 4 (a - l) 4 (a - 2) (2a - l)(a 2 - 3a + 3))u n 

— n 3 (n — a) 3 (3n - 2a) (3n + 1 - 2a)(3n - 1 - 2a) (39n 4 - 13(5a - 12)n 3 
+ 3(15a 2 - 65a + 78)n 2 - 3(5a 3 - 30a 2 + 65a - 52)n 

+ (2a 4 - 15a 3 + 45a 2 - 65a + 39))w n _i = 0, (10) 
uq = 1, u\ = (a — l)(a — 2)(a 2 — 3a + 3), 



2a 4 - 15a 3 + 45a 2 - 65a + 39 lim — = 6Z 4 (a), 

W = 0, Vl = — o , n ^°° U n 



and Theorem 2 in [22] yields the following 5-fold integral: 

(-l) n r(3n + 2-2a) 




- x x ) n n - =2 x]~ a (l - Xj ) n dxi • • • dx 5 



J (xi(l - (1 - (1 - (1 - x 2 )x 3 )x4)x 5 ) + (1 - x 1 )x 2 x 3 x 4: x 5 ) 
[o,i] 5 



(11) 
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4. Other recursions. Polynomial recursions for approximations to the numbers 

00 f_1\n-l 00 1 

z-(a) = y ^ — , z 2 ( a ) = y - ^ 

may be constructed by means of simpler (not well-poised) hypergeometric series. 
Namely, taking 

r n = u n Z-{a) -v n = (-ir Yi-l)'- 1 -, (t ~ 1)(t ~ 2) /" ( /~ n) r 

1V 1 K 1 ' (t- a)(t-a + l) ■■■(t-a + n) 



t=i 



1 x n ~ a (l - x) n , 
< — ; s — n — clx 

! Q (i + x) n+1 

we obtain the second-order recursion 

(n + l)(n + 1 — a)(2n — a)u n +\ — (2n + 1 — a)(6n 2 — 6(a — l)n + a(2a — S))u n 
+ n(n - a)(2(n + 1) - a)u„_i = 0, (12) 

uq = 1, u\ = — 2a + 3, vq = 0, i>i = -, lim — = Z^(a), 

a — 1 n^oo <U n 

while choosing 

n! • (t- l)(£-2)---(£-n) 



r n = u n Z 2 (a) -v n = (-l) n y 



7~t ((* - a)(t - a + 1) • • • (t - a + n))' 
y ' J J (l-xy)^ 1 

[0,1] 2 



we arrive at the recursion 



in 4 



(n + l) 2 (n + 1 — a) 2 (5n 2 — 6an + 2a 2 )u n+ i 

- (55n 6 - ll(16a - 15)n 5 + 2(117a 2 - 220a + 90)? 

- (160a 3 - 468a 2 + 388a - 85)n 3 + (56a 4 - 240a 3 + 316a 2 - 142a + 15)n 2 

- 2a(a - l)(4a 3 - 24a 2 + 32a - 9)n - 2a 2 (a - l) 2 (2a - 3))u„ 

- n 2 (n - a) 2 (5(n + l) 2 - 6a(n + 1) + 2a 2 )w n _i = 0, (13) 



1 1 _l q n 2a 2 -6a + 5 v n 
u — 1, ui = -2a + 3, fo = 0, vi = — ; —-z — , lim — = Z 2 (a). 

(a — \y n^oo u n 

The approach presented above allows to derive a higher-order polynomial recursions 
for simultaneous approximations to odd and even zeta values and their a-shifts (see 
[21])- 

5. Modular remarks. It is worth mentioning that Apery's recursions for ratio- 
nal approximations to £(2) and £(3) (that correspond to the case a = in (4) 
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or (13) and in (8)) have a very nice modular interpretation: the generating func- 
tion U(z) = J2^Lq u nZ n G %>[[z\] becomes a modular form after substituting a suit- 
able modular function z = z(t) (see [4, 7]). This phenomenon happens for several 
other Apery-like recursions as well (see [5, 16]); Zagier's technique in [16] allowed to 
construct new simple recursions for rational approximations to the numbers 

and Yl ^/T = 1 Z 2 (2 j ( Catalan ' s constant), 



where 



(=*)=!/ (mod 3) and = { 



for v even, 

v (mod 4) for v odd 



are quadratic characters. 

In spite of the complicated form of the recursions given in Sections 1-4 above, 
their solutions admit nice arithmetic properties if a is a rational number. For in- 
stance, the corresponding generating functions U (z) = Xl^Lo u n zU satisfy the prop- 
erty u(Az) G Z[[z]], where the integer A depends on a G Q (although a proof of the 
property in full generality for the recursions (4) and (10) is still beyond reach; see [10, 
14, 20] for particular results). Beukers' computations [6] show that it is hard to expect 
modular parametrizations except in the above mentioned cases of the recursions for 
rational approximations to £(2) and C(3): the (Zariski closures of the) Galois groups 
associated to the linear differential operators annihilating functions U(z) are richer 
than SL 2 - Beukers considers the recursion for rational approximations to ((4) (the 
case a = in (10)) and shows that the corresponding differential Galois group turns 
out to be O5, while the linear differential operator corresponding to the recursion for 
rational approximations to Catalan's constant (the case a = 1/2 in (4)) is reducible. 
We expect that this differential reducibility holds for all a ^ Z and that irreducible 
components of the corresponding differential operators are pullbacks of hypergeomet- 
ric differential operators. The latter fact is closely related to the general conjecture 
(due to Dwork, Bombieri, . . . ) on the structure of differential G-operators. While 
chances to be able to attack this general conjecture seem to be small at the moment, 
to us it appears to be a nice and quite realistic program to give a direct proof in the 
case of the above recursions. 
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